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ABSTRACT 

X-ray observations of synchrotron rims in supernova remnant (SNR) shocks show evidence of ef- 
ficient electron acceleration and strong magnetic field amplification (a factor of ~ 100 between the 
upstream and downstream medium). This amplification may be due to plasma instabilities driven 
by shock-accelerated particles, or cosmic rays (CRs), as they propagate ahead of the shocks. One 
candidate process is the cosmic ray current-driven (CRCD) instability (Bell 2004), caused by the elec- 
tric current of "unmagnetized" CRs (i.e., CRs whose Larmor radii are much larger than the length 
scale of the CRCD modes) propagating parallel to the upstream magnetic field. Particle-in-cell (PIC) 
simulations have shown that the back-reaction of the amplified field on CRs would limit the amplifi- 
cation factor of this instability to less than ~ 10 in galactic SNRs (not including the additional field 
compression at the shock). In this paper, we study the possibility of further amplification driven near 
shocks by "magnetized" CRs, whose Larmor radii are smaller than the length scale of the field that 
was previously amplified by the CRCD instability. We find that additional amplification can occur 
due to a new instability, driven by the CR current perpendicular to the field, which we term the 
perpendicular current-driven instability (PCDI). We derive the growth rate of this instability, and, 
using PIC simulations, study its non-linear evolution. We show that the maximum amplification of 
PCDI is determined by the disruption of CR current, which happens when CR Larmor radii in the 
amplified field become comparable to the length scale of the instability. We find that, in regions close 
to the shock, PCDI grows on scales smaller than the scales of the CRCD instability, and, therefore, 
it results in larger amplification of the field (amplification factor up to ~ 45). One possible observa- 
tional signature of PCDI is the characteristic dependence of the amplified field on the shock velocity, 
B 2 oc v 2 h , which contrasts with the one corresponding to the CRCD instability acting alone, B 2 oc v^ h . 
Our results strengthen the idea of CRs driving a significant part of the magnetic field amplification 
observed in SNR shocks. 

Subject headings: ISM: magnetic filed — cosmic rays — supernova remnants — jets and outflows 



1. INTRODUCTION 

Cosmic rays (CRs) up to ~ 10 15 eV are most likely ac- 
celerated in galactic supernova remnant (SNR) shocks 
via the diffusive shock acceleration (DSA) mechanism 
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ord fe Ostriker||1978p . This mechanism is 
be highly non-linear, due to the expected coupling be- 
tween the evolution of the thermal plasma, the CRs, and 
the magnetic turbulence in the shock vicinity. Indeed, 
there is now observational evidence suggesting that elec- 
trons are being accelerated in the non-relativistic shocks 
of SNRs, and that significant magnetic turb ulence is be- 
ing produced as part o f the process (see, e.g. |Ballet|2006 
Uchiyama et al. 2007). The evidence comes from X-ray 
observations ot SlNRs that show the existence of thin, 
non-thermal rims, which are interpreted as synchrotron 
emission by TeV electrons accelerated at the shocks. The 
rapid variability and thinness of the rims (which depend 
on the synchrotron cooling time of the electrons) have 
allowed to estimate the strength of the field, suggesting 
downstream amplitudes ~ 100 times larger than typi- 
cally expected in the ISM of the Galaxy. This implies 
a significant amplification even accounting for the com- 
pression of the field at the shock, which would contribute 
a factor of ~ 4 to the growth of the field. 
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The origin of this strong amplification is an open ques- 
tion. On the one han d, using MHD shock simulations, 



Giacalone et al. ( 2007 I showed that the presence of den- 
sity and magnetic fluctuations in the upstream medium 
of non-relativistic shocks may introduce vorticity into 
the flow and p roduce field amplification at the shock 



Sironi & Goodman||2007|). This mecha- 



itself (see, also, 
nism would produce a magnetic enhancement larger than 
expected from the simple shock compression of the up- 
stream field. On the other hand, it has been proposed 
that the magnetic field may be amplified by plasma in- 
stabilities driven by the CRs themselves, as t hey propa - 
gate through the upstream medium of shocks ( |Bell|2004| . 
This idea is interesting because, besides helping to ex- 
plain the field amplification, the implied magnetic tur- 
bulence would increase the CR confinement to the shock 
vicinity, enhancing the efficiency of DSA. Due to this, 
understanding the extent to which magnetic fields are 
amplified by CRs, as well as characterizing the corre- 
sponding magnetic turbulence, is essential for any realis- 
tic model of DSA. 

The cos mic ray curre nt-driven (CRCD) instability pro- 



posed by Bell (2004) has gathered considerable atten- 
tion recently. The CRCD instability consists of circu- 
larly polarized Alfven-type waves, driven by the elec- 
tric current of CRs propagating along the magnetic field 
lines. It requires the CR Larmor radii, i?L iCr , to be 
much larger than the wavelength of the fastest grow- 
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ing mode, X C rcd = cB /J cr , where B , J cr , and c 
are the magnitude of the initial magnetic field, the CR 
current, and the speed of light, respectively. This fea- 
ture makes this instability different from the previously 
proposed resonant instability, in whi ch the waves grow 
at wa velengths comparable to Rl,ct ( |Kulsrud fc Pearce 



1969 1. The growth rate of the CRCD instability is given 



2004), where p is the 



by icrcd = J cr (ir/pc 2 ) 1/2 peU| 

mass density of the background plasma. For the typical 
conditions of CR-modified SNR shocks (where the en- 
ergy densities of CRs, u cri downstream thermal plasma, 
u t h, and magnetic field, ub, satisfy u t h > u cr 3> ub), this 
rate is expected to be significantly larger than that of the 
resonant instability. 

The non-linear evolution of t he CRCD instability has 
be en stu died using both MHD ( |Bell[2 004; ZirakashvUTet] 
aT1|2008D and particle-in-cell (PIC) ([N icmicc et al.l|2008| 
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ulations. rioth kinds of studies showed that the CRCD 
instability has the potential of growing to very non-linear 
amplitudes as long as the CR current is kept constant, 
i.e., if the back-reactio n on the CR trajectories is ne- 
glected. In particular, jRiquelme fc Spitkovskyj (2009) 
found that, at constant CR current, the field growth 
stops when the Alfven velocity of the background plasma, 
va, gets close to the CR drift velocity, 

v d,cr- However, 

PIC studies show that if the CR back-reaction is con- 
sidered, the field can also saturate due to the trapping 
of the CRs in the amplified field. When Rl,ct becomes 
comparable to the dominant length scale of the magnetic 
turbulence, the CRs get strongly deflected, which signif- 
icantly reduces their current and quenches the magnetic 
growth. This is the most likely saturation mechanism in 
the case of SNRs, which implies a maximum amplification 
factor of ~ 10 in the upstream medium of shocks, and 
confirms the CRCD instability as a viable mechanism 



for magnetic turbu lence generation in SNRs (Riquelme 
& Spitkovsky||2009|) . When the field compression at the 



shock is considered, this amplification may increase by an 
extra factor of ~ 4. However, since the upstream amplifi- 
cation factor of 10 constitutes an upper limit, the CRCD 
instability alone would not be enough to account for the 
factor of ~ 100 inferred from the X-ray observations. 
In this paper we study the possibility of magnetic ampli- 
fication beyond the saturation of the CRCD instability. 
The situation we explore is one where CRs propagate 
through a plasma where previous CRCD magnetic tur- 
bulence has already been produced on scales larger than 
the typical CR Larmor radius, Rl,cv Thus, our analysis 
would be applicable to the lower energy CRs, which are 
more confined to the shock vicinity. This is motivated 
by the fact that the CRCD instability is probably driven 
first by the highest energy CRs, in regions far upstream 
from the shock. This instability saturates when the Lar- 
mor radius of these high energy particles is about the 
length scale of the pre-amplified magnetic fluctuations, 
Ao- Thus, as the shock approaches, a large fraction of 
the CR energy will be carried by "magnetized" CRs, i.e., 
CRs with Rl,ct smaller than Ao- Then, near the shock, 
CRs will find regions of pre-amplified, quasi-transverse 
field, which can be considered homogeneous on scales 

~ Rh.cr- 

If Rl,ct ;$ Aq, then CRs will not easily diffuse through 



the regions of pre-amplified field. Instead, they will pro- 
trude into these regions only by a distance of ~ Rl,ct- 
This situation will produce a CR current, J cr , perpen- 
dicular to the initial, pre-amplified field, B 0l due to the 
coherent deflection experienced by CRs in the "homoge- 
neous" (on scales of ~ Rl,ct) magnetic field. We propose 
that, under these conditions, an extra magnetic amplifi- 
cation is possible due to a new instability: the perpen- 
dicular current- driven instability (PCDI), which consists 
of purely growing, compressional waves that arise when 
J cr is perpendicular to Bq. 

As we show below, the PCDI would produce larger mag- 
netic amplifications compared to the CRCD instability 
acting alone. In addition, this instability can amplify 
magnetic fluctuations on scales comparable to the Lar- 
mor radii of the lowest energy CRs, which would im- 
prove their diffusion and increase the efficiency of their 
acceleration at the shock. If this did not happen, the 
large-scale, transverse fields generated by CRCD insta- 
bility would hamper the low energy CR diffusion, which 
would probably shut off the acceleration of the CRs that 
cause the CRCD amplification in the first place. 
In Sp] we explain the physics of the PCDI. First, we 
showhow the penetration of CRs into the regions of pre- 
amplified field can give rise to a current, J cr , perpen- 
dicular to the pre-amplified field, B Ql and why we expect 
that situation to happen in the upstream medium of SNR 
shocks. Then, we show how the presence of a Jcr perpen- 
dicular to Bq can produce the PCDI, and derive its dis- 



persion relation. In \ 
of the PCDI using P 



f}3l we show the results of our study 
simulations. First, we model the 
non-linear evolution of the instability assuming a con- 
stant CR current, i.e., ignoring the back-reaction on the 
CRs, and provide a simple analytical model for the non- 
linear behavior of the PCDI modes. Second, we show 
the results of a series of simulations that model the non- 
linear evolution of the PCDI including the full CR dy- 
namics. The main focus is to show the role of the CR 
back-reaction on the final saturation of the instability. 
In p] we apply our results to the case of SNR shocks. 
Finally, in f|5]we summarize our results and present our 
conclusions. 

2. PHYSICS OF PCDI 
The perpendicular current-driven instability (PCDI) is 
caused by a CR current, J c r , pe rpendicular to the ambi- 



ent magnetic field, B . In { 2.1 we show why we expect 
this J cr _L Bo to actu ally occur in the precursor of SNR 
shocks. Then, in §2.2| we explain the physics of the PCDI 
growth, and derive its dispersion relation. 

2.1. Perpendicular Current in SNR Shock Precursors 
In this section, we show that a significant CR mean ve- 
locity perpendicular to an initial field, _B n , can arise in 
the upstream medium of SNR shocks, if there is a sig- 
nificant anti-correlation between the CR number density, 
n cri and the strength of the magnetic field. 
If CRs propagate through a smooth medium, as may be 
the case shortly after the shock formation, there is no rea- 
son to think that a significant anti-correlation between 

n cr and \Bq\ can happen. However, if some magnetic 



Perpendicular Current-Driven Instability 



3 



downstream medium 



•. CRs 

... ■ ■ . 



. shock c ; .• ■ 



y ■ 



pre-amplified field, B ( 




R.egian"A" 



ptt- amplified field, B ® 



* * ► I* ■ r I ! 

t 1 * 1 l 1 *■ * 

■* I I ■ i I ■ 

M i a a ■ « l a < * * a • • 

: . I« J» :•■ «• Kj 

i r * ■ . . -f 

. crs , ■■ ■ • ./T7 . . 

:.. : - .< i . ■.. ■.. • .< 



&; ■ .'-> : 

y :i " : x • 
I* 




Region"B" 



Fig. 1. — A diagram with the global picture of formation of zones 
with CR current perpendicular to the pre-amplified magnetic field, 
Bq, m the upstream medium of SNR shocks. The upper panel 
shows the downstream and upstream regions of a generic shock. 
Regions of pre-amplified field in the upstream are represented by 
grey, doughnut-like shapes. Black dots represent "magnetized" 
CRs, i.e., CRs whose Larmor radii are smaller than the size of the 
regions of pre-amplified field. The lower panel shows a zoomed 
and rotated (notice the orientation of coordinate axes in both pan- 
els) image of the small, solid-line rectangle shown in the upper 
panel. This image is divided into region "A" (with low initial mag- 
netic field), and region "B" (with the pre-amplified magnetic field 
So)- The CRs, coming from region "A", can penetrate into region 
"B" by a distance ~ 2Rl o, which defines the interface into which 
CRs can protrude. As indicated in the figure, the coherent deflec- 
tion of the CRs in this region produces a net mean velocity along 
Z, of magnitude ~ c/2. Finally, the small, dashed-line rectangle 
depicted in th e upper panel represents the type of regions that we 
study in j]3.3|by using two-dimensional PIC simulations. 



turbulence has already been amplified, then this anti- 
correlation can exist. For instance, numerical MHD and 
PIC simulations of the CRCD i nstability (see Fig. 4 of 
|BellN2004| or Figs. 10 and 12 of |Riquelme fe Spitkovsky] 
( |2009| ) show how this instability produces turbulence 
characterized by typical plasma density contrast of ~ 10, 
with frozen-in magnetic field concentrated in the high 
densit y regions. As discussed in Riquclmc & Spitkovsky 
( 2009 ) , this kind of turbulence can be produced far up- 
stream from the sh ock by high-energy CRs "escaping" 
at ~ c/2 (see also |Zirakashvili et al. 20081. This im- 
plies that lower energy CRs, which are more confined 
to the shock, will encounter the upstream medium with 
previously amplified magnetic turbulence. This is the 
situation depicted in the upper panel of Figure [TJ which 
shows the regions in the upstream of a shock contain- 
ing a pre-amplified magnetic field, Bq (diagrammatically 
represented by doughnut-like shapes in the figure; note 
that, in reality, these "doughnuts" are not necessarily in 
one plane) surrounded by CRs (indicated by black dots). 



Since the saturation of the CRCD instability happens 
when the Larmor radius of the highest energy CRs is 
close to the length scale of the pre-amplified field, Ao, 
then most CRs near the shock are expected to have Lar- 
mor radii, Rl,ct, that satisfy Rl.ct < ^a- We will refer 
to CRs satisfying this condition as "magnetized" CRs. 
In this scenario, CRs will not be able to diffuse through 
the regions of pre-amplified field. Instead, they will only 
be able to penetrate into these regions by a distance of 
~ 2Rl ct . The lower panel in Figure [TJdepicts this situ- 
ation. It shows a zoomed image of the interface between 
a region of low and high upstream magnetic field (rect- 
angular box in the upper panel). This plot is divided into 
two regions. Region "A" (left), which contains the CRs 
and where there is almost no initial magnetic field, and 
region "B" (right), where there is a pre-amplified mag- 
netic field, Bq, which has a length scale, Ao, larger than 
the CR Larmor radii, Rl. cr . Since Rl.ct ^ Ao, CRs pro- 
truding into region B (by a distance ~ 2i?r J . cr ) will not 
experience random scatterings due to fluctuations in the 
field. Instead, they will just complete a fraction (typi- 
cally ~ 1/2) of a Larmor gyration before being scattered 
back into region A. The fact that CRs complete only 
a fraction of a Larmor cycle in region B implies that, 
while in this region, they will tend to move parallel to 
Vn cr x Bq (i.e., along z in Figure [lj. Thus, if CRs have 

positive electric charge, they will produce a current, J cr 
perpendicular to the pre-amplified field 1 . Since CRs are 
relativistic, their mean velocity perpendicular to Bq will 
be a fraction (of the order of ~ 1/2) of the speed of light, 
larger than the typical drift velocities along the field that 
drive the CRCD instability. This suggests that it is pos- 
sible to find upstream regions near SNR shocks where J cr 
can be quasi-perpendicular to Bq. The formation of the 
per pend icular mean CR velocity at ~ c/2 is confirmed 
in §3. 3| using two-dimensional PIC simulations. 



2.2. PCDI Dispersion Relation 

The presence of the CR current, J cr , perpendicular to 
the initial field, Bq, induces a return current, J re t, al- 
most equal to — J cr in the background 2 . Because of this, 
a constant force, — J cr x Bq, will push the plasma ho- 
mogeneously. Thus, if plasma density fluctuations, Sp, 
exist, this — J cr x Bq force will accelerate the low density 
regions more than the ones with higher density. As we 
will see below, this differential plasma acceleration in the 
direction perpendicular to Bq can amplify the transverse 

component of the magnetic field, SB, and enhance the 
density fluctuations in the plasma, producing the expo- 

1 Since negative CRs would move in the opposite direction, their 
negative charge would finally make their current point parallel to 
the one of positive CRs, implying an even larger J cr . 

It is interesting to note that, even though charged particles 
can not easily move across magnetic field lines, in the case when 

J cr X Bo, J re t is still almost equal to —J cr - In this case the 
compensating current is due to the polarization drift of particles, 
v p = (mc 2 / eB^dE/ dt, where m and e are the particle mass and 
charge, and E is the electric field. It is possible to show from the 
Ampercs's law that, if t>A,o/ c 1, the return current produced by 
this drift (acting mainly on the ions) almost completely compen- 
sates the CR current (see also Lyutikov 2009). 
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nential growth of the PCDI modes. We note that this 
instability starts from a non-equilibrium state, where 
the whole plasma is being accelerated in the direction 
of — J cr x Bq. 

We derive the dispersion relation of the PCDI in the 
MHD limit. We assume the existence of an externally 
imposed CR current, J cr _L Bq. The way this CR current 
is produced and the back-reaction of the amplified field, 
B (= B + SB), on the CR trajectories is not modeled in 
this derivation. The current carried by the background 
plasma is cV x B/Att — J cr , which includes the return 
current induced in the background. So the MHD equa- 
tions will only be modified by an extra term in the force 
per unit volume on the plasma: 



dp 

at 



+ V • (pv) = 



(1) 



,dv 



p(^ + (w-V)u) + Vp- — (VxB)xB+- xB = (2) 
at iir c 



dB 
~dt 



Vx(»xB) = ! 



(3) 



where v is the plasma fluid velocity and p is its pressure. 
Although the CRs can also include accelerated electrons, 
which would reduce the CR net charge, in general CRs 
may carry a charge density, p c>cr , that must be com- 
pensated by the background plasma. This extra plasma 
charge should introduce an electrostatic term in Equa- 
tion p| due to plasma charge density fluctuations. It is 
possible to show that this term is second order in per- 
turbation theory, so it can be neglected in our linear 
analysis. We refer the reader to Appendix [B] for further 
details. 

If each physical quantity £ is separated into an initial, 
homogeneous component, £o, plus a small perturbation, 
< £o), (so B = B +SB, v = vq+5v, andp = p +Sp), 
then the MHD momentum equation will have a zeroth 
order part: 

dv 

which represents homogeneous acceleration of the 
plasma. The first order quantities read: 



Jcr x Bq — 0, 



(4) 



dSp 



+ p V ■ (Sv) = 



(5) 



Per 



dSv 
dt 



c 2 s VSp 
1 



SpJcr 

Po c 



X Bq 



47T 



(V X SB) X Bq 



x6B = (6) 



dSB 
~0T 



- V x (Sv x B a ) = 0, 



(7) 



with c s being the sound speed in the plasma and where 
we have chosen a reference frame where \vo(t = 0)| = 0. 
The physical meaning of the different force terms in 
Equation ^ can be understood with the help of Figure 

[ij Let perturbations SB and 5p be such that the mag- 
netic field lines and plasma density are represented by 



the solid and dotted lines of Figure [2] respectively. Then, 
since the force J ret x Bq is the same everywhere, regions 
of low density will accelerate more than regions of high 
density, making the overdense regions lag behind with 
respect to the rest of the plasma. This is the origin of 
the force Spj PqJ ct x Bq in Equation (roj. As shown by the 
solid black arrows in Figure [2] this force tends to stretch 
the lines of magnetic field and, therefore, produces the 
amplification of SB. This amplification can happen ex- 
ponentially, because the growth of SB increases the force 
~J cr x SB, which, as depicted in Figure [i] enhances the 
density contrast even more. In a sense, the PCDI is 
similar to Parker instability, since the magnetic growth 
implies the existence of regions of the plasma that accel- 
erate more than others in the direction perpendicular to 
Bq. However, whereas in the case of Parker instability 
the amplification is driven by a homogeneous gravita- 
tional acceleration, the amplification of PCDI is caused 
by a homogeneous force, given by — J cr x Bq. 
Performing the Fourier transform of Eqs. and 
|7| , and assuming that the Fourier modes are stationary, 

with fe 1 1 -Bq , it is straightforward to obtain the dispersion 
relation: 



7 



+ c 2 k 2 



Po(l 2 + c 2 s k 2 ) ' 



(8) 



where 7 is the growth rate corresponding to a given k. 
This dispersion relation can also be obtained by taking 
the limit k\\B _L J „ r of the m ore general MHD result 



shown in Eq. (4) of Bell ( 2005 1 . A derivation using the 
multifluid approach can be found in Appendix |A]) . From 
Equation ([8]) we obtain that the fastest growing mode has 
the wave number, 



^pcdi — 



AirJ cr 
c 2 B Q I 



' J Xo/ C s _C 



(9) 



and the growth rate, 



Ipcdi - 2J cr (^c 2 ) 1 / 2 i ^ o/e - / , (10) 



where va,o(= Bq/^/A-kpq) is the initial Alfven veloc- 
ity in the plasma. Thus, if the background plasma 
is initially cold (c s — > 0), then kpcDi — * °o and 
Ipcdi ~ 2J cr (ir / pc 2 ) 1 / 2 . However, if we assume that 
c s is close to the typical turbulent velocity in the 
plasma, then it is possible to show from Eqns. J5|7] 
that the sound speed due to PCDI turbulence would be 
~ va,oSB I Bq. Thus, the thermal effects would make 
kpcDi ~ 2nJ cr /cSB and would keep the same growth 
rate jpcdi ~ 2J cr (-K / pc 2 ) 1 / 2 as long as c s < va,o- 
As SB/Bq — > 1, the expressions for the wave number 
of fastest growth, kpcDi, and its corresponding growth 
rate, jpcdi, look the same as the ones for CRCD insta- 
bility: kcR CD = 2TrJ cr /cB and JcrCd = Jcr(n / pc 2 ) 1/2 
(Bell 2004). However, the CR currents that drive the 



PCDI and CRCD instabilities are the components of J cr 
perpendicular and parallel to -Bo, respectively. Thus, if 
J cr is quasi-perpendicular to Bq (as may be the case in 
regions close to the shock) , the PCDI would grow faster 
than the CRCD instability. 
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Fig. 2. — A diagram illustrating the physics of PCDI. If an ex- 
ternal current J cr is applied perpendicular to the plane of the fig- 
ure and also perpendicular to the initial field, Bo, the background 
plasma will develop a return current. J re t ~ — Jcr pointing out 
of the plane of the figure. This will produce a homogeneous force 
— Jcr x Bo accelerating the plasma, as shown in the diagram. The 
PCDI begins from this non-equilibrium state. Let us consider the 
magnetic field lines and the plasma density are represented by the 
solid and dotted lines of the figure, respectively. Then, since the 
force —J cr x Bo is the same everywhere, regions of low density will 
accelerate more than regions of high density, making the overdense 
regions lag behind with respect to the rest of the plasma. This is 
the origin of the force Sp/poJcr x Bo of Eq. jfi]l (shown by the solid, 
black arrows), which tends to stretch the lines of magnetic field and 
amplify SB. The other forces from the same equation: — J cr x SB, 
— c^\7Sp, and (V X SB) X Bo/4ir, are represented by the solid grey 
arrow, the dotted grey arrow, and the dotted black arrow, respec- 
tively. For unstable modes, this amplification happens because the 
growth of SB increases the force — J cr x SB, which enhances the 
density fluctuations. Thus, the growth of SB increases the force 
that produces it, resulting in an exponential growth of SB and Sp. 

3. PIC SIMULATIONS 

In this section we use two-dim ensi onal PIC simulations 
to study the PCDI. First, in §3.2[ we explore the case 
of a constant, externally imposed CR current, J cr , p er- 



pendicular to the ambient field, B . Then, in £3.3 we 
study the full dynamical evolution of the CRs, allowing 
a non-constant J cr . We pay special attention to the way 
a current J cr perpendicular to Bg can set in, and to the 
effect of the CR back-reaction on the saturation of the in- 
stability. In each section we specify the numerical set up 
of the corresponding simulations, listing their common 
features in 93.11 



3.1. Common Simulation Setup 

We use the electromagnetic PIC code TRISTAN-MP 
(Buncman 1993; Spitkovsky 2005). Our simulations use 
periodic, two-dimensional boxes whose normal vector is 
chosen to point along z. The boxes contain a charge- 
and current-neutral plasma, consisting of cold ions and 
electrons (with initial thermal velocities ~ 10~ 4 c), which 
constitute the background plasma, and a population of 
relativistic particles or CRs that may either be composed 
of ion-CRs (positive CRs) only, or a mixture of both 
ion- and electrom-CRs. The boxes also contain the initial 
magnetic field, Bq, in the plane of the simulation, point- 
ing along the y axis. Particles are spread within the 



box in such a way that the local charge density is zero 
everywhere. Since typically in a SNR the CR density 
much smaller than the background ion density, rij, 
we increase CR statistics by modifying the charge of the 
plasma and CR macroparticles (for instance, by making 
Qi,cr ~ Qt n i,crl n ii where Qi and Qi. cr are the charge 
of the ion and ion-CR macroparticles) . The mass of the 
macroparticles is modified accordingly in order to keep 
their charge to mass ratios unchanged. 

3.2. Constant CR current 

In this section we focus on the case of a constant J cr 
pointing along z (out of the plane of the simulation). 
Thus, the back-reaction of the field on the CR trajec- 
tories is explicitly suppressed in these simulations. We 
show the results of four runs with different n cr /rii ratios 
and CR charge densities, so that we can test the disper- 
sion relation (Eq. [8]) and investigate the effect of hav- 
ing different CR net charge. In two of our simulations 
the CRs have no net charge (they are made of equal 
amount of ion- and electron-CRs), while in the other 
two runs CRs are only positively charged. In the case 
of zero CR net charge, the CR current is produced by 
two counter-streaming CR beams of equal density but 
opposite charge. For both the zero and non-zero CR 
charge cases we test n cr /ni = 0.01 and 0.04, where n cr 
is the total density of CRs (which in the case of zero CR 
charge includes the positive and negative CRs) 3 . In both 
types of simulations, particles are randomly distributed 
in the box, but placed such that, initially, there is lo- 
cal charge neutrality. Also, the initial local current is 
made zero by giving electrons a small initial velocity so 
that they carry a compensating current —J cr . The other 
physical and numerical parameters of the simulations are: 
i>a,o/c — 1/40, mass ratio mi/m e = 10, CR current ve- 
locity Vd.cr = 0.6c, speed of light c = 0.45A/At, the 
skin depth of electrons c/uj p . e = 3A, and the number of 
macroparticles per cell per species is 16. Here, A is the 
grid cell size and At is the time step. 
First, we tested the dominant length scale for the PCDI 
amplification. Figure [3] shows the space distribution of 
the field and background plasma density for the run with 
ricr/ni — 0.01 and neutral CR beam at three different 
times [typcDi = 17, 21, and 35, where jpcdi is the theo- 
retical growth rate of the simulation) . When SB / B < 1 
(tjPCD = 17, in panel a), the dominant wavelength is 
smaller than Xpcdi (which is chosen to be the wave- 
length of the fastest growing mode when SB /Bo = 1). 
When 5B/B « 1 (tjpcD = 21, in panel b), the length 
scale of the fluctuations becomes very close to Xpcdi, 
confirming our analytical result. As the field becomes 
non- linear, it adopts the shape of loops, which tend to 
merge and grow in size. We see that in the very non- 
linear regime (tjpcD = 35) the dominant scale of the 
loops grows by a factor comparable to the magnitude of 



3 In the zero CR charge case, the charge and number densities 
of macroparticles 

-(Qc/2)n C r/rii, and 7V;. cr = N e .cr = Ni 



N e (where Nj rep- 



resents the number of macroparticle per cell for species "j"). In 
the non-zero case, Q cr = —Qen-cr/ni- Then, in order to have 
local charge neutrality we need Qi = —Q e (l — n cr /ni), and 
N cr = Ni = N e (this way there is a small excess of electron charge 
in the background plasma, which compensates the charge of the 
CRs). 
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the amplification, which is consistent with our predic- 
tions. 

Our simulati ons also confirm the analytic growth rate, 
given in Eq. (10 1. Figure El shows the evolution of the 
energy in the three magnetic field components in our 
simulations. The left and right panels represent the cases 
with n cr /rii — 0.01 and 0.04, respectively. Solid and dot- 
ted lines show the cases with zero and finite CR charge. 
The colors black, green, and red represent the magnetic 
energy of the x, y, and z field components. The time is 
measured in units of ^pcpj, where 'fpcDi is the PCDI 
theoretical growth rate for n cr /rii — 0.04. In the linear 
regime, the PIC results are consistent with our theoret- 
ical estimates for the PCDI growth rates. Also, we see 
that, as long as SB -C Bq, the growth rate is indepen- 
dent of the net CR charge. This situation changes in the 
non- linear regime. While in the zero CR charge cases 
the magnetic field keeps growing with about the same 
growth rate as in the linear regime, the non-zero CR 
charge simulations saturate at 8B/B ~ 1. 
The non-linear behavior of PCDI can be understood 
as follows. As the amplified field becomes non-linear, 
it adopts the shape of magnetic loops that surround 
background plasma holes on scales of ~ Xpcdi (as 
shown in Figure pi). The evolution of magnetic loops 
in the presence of a constant, perpendicular CR cur- 
rent w as already considered by |Milosavljevic fc Nakar 
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(20061, whose argument we revise here. If J cr is con- 
stant, the background plasma will create a return cur- 
rent, J ret rs — J cr . Thus, this return current will pro- 
duce a force w —J cr x B that tends to expand the loops. 
In the MHD regime, as a loop expands, the magnetic 
flux freezes in the fluid, keeping the ratio B/pr con- 
stant (where r is the radius of the loop). Thus, if no 
extra forces are considered, we obtain that the acceler- 
ation of the loop is d 2 r/dt 2 ~ J cr Bor/\pcDiPoc, which 



is solved by r(t) ~ \pcdi& 1fcl , where Xpcdi is equal 
to 2ir /kpcDi in the limit 5B = Bq. Thus, since flux 
freezing keeps B/pr constant, the magnetic field should 
also grow exponentially at approximately the same rate 
as r{t). However, we see that, when CRs are made of 
only positively charged particles, this non- linear ampli- 
fication gets suppressed. In that case, the expansion of 
the negatively charged loops would produce charge sep- 
aration in the plasma. Due to this, a monopolar electric 
field will form in the loops, which would suppress their 
expansion. Although this force plays an important role 
in the cases analyzed in this section, in cases where J cr is 
caused by CRs protruding into regions of pre-amplified 
field (like in the lower plot of Figure [TJ , they would not 
stop the growth. The reason has to do with the particu- 
lar way the charges are distributed within each loop. We 
will analyze in detail this point in Appendix |B] 
In summary, in this section we have used two- 
dimensional PIC simulations to show that, if the CR 
current is kept constant, the PCDI grows at the wave- 
lengths and rates predicted by the linear theory. In the 
non-linear regime, the exponential growth can continue 
with about the same rate as in the linear regime, and 
with a length scale proportional to SB. However, if CRs 
are mainly positively charged particles, this growth may 
be affected by electrostatic fields due to charge separa- 
tion in the plasma (see Appendix [B| . 




Fig. 3. — The space distribution of the magnitude of the gener- 
ated field and the background plasma density at typcDl = 17, 21, 
and 35 (represented by panels a), b), and c), respectively), for 
the run with n cr /rii = 0.01 and neutral CR beam. The initial 

magnetic field, Bo, and the CR current, J cr , point along y and z, 
respectively. The distance is normalized in terms of the theoretical 
XpCDl a t <5B /Bo = 1. We see that in the linear regime (panel 
a) the size of the fluctuations, A, is smaller than Xpcdi, and that 
it evolves into A ~ (5B / Bo)Xpcdi as the amplification becomes 
non-linear (panels b) and c). When SB/Bo 3> 1 (panel c), the 
PCDI modes acquire a loop-like shape. 
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Fig. 4. — The evolution of the magnetic energy in the three carte- 
sian axes (black, green, and red lines correspon d to the x, y, and 
z axes, respectively) for the four simulations in j ]3.2| The left and 
right panels show the cases with n cr /rii = 0.01 and 0.04, respec- 
tively. The solid lines show the situations when the background 
plasma (and the CRs) have a zero net charge, while the dotted lines 
represent the cases when the background plasma has a negative net 
charge compensated by the positive CRs. The time is measured in 
units of 'ypQjjj , where ypcDI ls f ne theoretical growth rate of the 
instability for the case n cr /nj = 0.04. 

3.3. CR back-reaction 

In this section , we in clude the full dynamics of CRs our 
the study. In §3.3. 1[ we show how CRs protruding into 
regions of pre-amplified field can produce a J c r _L B q, 



and drive the growth of PCDI modes. Then, in { 3.3.2 we 
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study the effect of the CR back-reaction on the saturation 
of the instability. 

3.3.1. Growth of PCDI 

In our simulations we want to reproduce the situation 
pictured in Figure [T] where regions with a pre-amplified 
field (represented by grey "doughnut-like" shapes) are 
surrounded by a population of magnetized CRs (i.e., CRs 
with Larmor radii smaller than the length scales of these 
regions). In particular, we want to model small sections 
of these regions, like the one marked with a dotted-line 
box in the upper panel of Figure [T] The idea is to study 
the growth of the PCDI due to CRs protruding into 
the regions of pre-amplified field region. Since we use 
a monoenergetic population of CRs, these particles are 
expected to penetrate by a distance comparable to two 
times their Larmor radius, i?L jCr , which defines a region 
or interface where the field grows. The other parts of 
the plasma, which are not reached by these CRs, should 
experience compression due to the acceleration of the 
interface. As CRs turn in the field of the magnetized 
region, the acquire the mean velocity in the z direction 
(out of the box), and the associated perpendicular cur- 
rent drives the rippling of the field characteristic of the 
PCDI. 

The numerical setup we use is similar to t he o ne o f the 
constant CR current case described in §3.1| and |3.2| The 
only difference is that, in this case, CRs are positively 
charged and monoenergetic particles, and are isotrop- 
ically injected within a narrow strip along the y axis. 
This strip is located in the left part of the simulation 
box and its width is l/40th of the width of the box. 
Thus, whereas the strip represents the region with high 
CR density and low initial magnetic field, the rest of the 
simulation box would represent the inner part of the re- 
gion of pre-amplified field. Extra electrons are injected 
along with the CRs so that the net injected charge is zero. 
The injection criterion is such that the mean density of 
CRs in the strip, n crj o, is kept constant 4 . 
We ran a series of simulations that share the follow- 
ing parameters: mi/m e = 10, n cr a /n i a = 0.04 (where 
riifi is the initial background ion density in the injection 
strip), T cr = 5, c = 0.45A/A£, c/w Pje = 6.4, and use 
8 macroparticle per cell per species. The difference be- 
tween the runs is given by their initial Alfven velocity, 
va,o/c, which ranges from va,o/c — 1/160 to 1/40. The 
evolution of the relevant physical quantities for a sim- 
ulation with va,o/c — 1/80, which we will call "5*", is 
depicted in Figures [5] M and [71 These figures correspond 
to t'jpcDi — 32,40 and 44, where ^pcdi is the growth 
rate of the PCDI considering n cr = n crj o and Vd.cr = c/2. 
In panels a) of these Figures, we depict the space struc- 
ture of the generated field. As CRs penetrate into the 
regions of zero initial CR density, we see the formation 
of magnetic loop s, sim ilar to the ones shown in the case 
of constant J cr ([3.2 1 5 . In panels b) we show the density 
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4 Similarly to the case with constant CR current, the charge of 
the macroparticles is modified such that it is possible to have a large 
number of CR macroparticles per cell. So, Q cr = Qin cr fi/ni,o, 
where n, q is the physical density of ions in the strip at t = 0. 

5 Since the simulation box is periodic, CR propagate both to 
the left and to the right of the injection strip. Thus, two nearly 
symmetric regions of field amplification form on both ends of the 
box. 
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Fig. 5. — Distribution of different particle and field quantities 
for simulation S* at t7pcu/ = 32, where ■ypcDl is the growth 
rate of the PCDI calculated using n cr o and c/2 as the number 
density and mean velocity of CRs. The distance is shown in units 
of R° L , which is the Larmor radius of the CRs in the initial field 

So (= Boy). The simulation box is periodic, a) Magnitude of the 
generated field; b) Density of different species averaged over y as 
a function of the x coordinate. The red and green lines represent 
the background ions and electrons, whose densities are normalized 
by «i,o> an d the black and blue lines represent the densities of 
CRs and injected electrons (to compensate CR charge), which are 
normalized in terms of n cr o; c) Mean velocities of CRs along x 
(black line) and z (blue line), normalized by c; d) Magnetic field 
amplification as a function of x, normalized in terms of Bo- The 
red and blue lines show the maximum and average amplifications, 
respectively. Black line shows the field magnitude corresponding to 
local energy equipartition with the CRs. Dashed lines in panels b), 
c) and d) show the results of a simulation analogous to S* , but with 
the y dimension of the box significantly reduced (quasi-lD). This 
simulation shows only the effect of compressional amplification. 



of the different species as a function of the x coordi- 
nate. It shows how CRs can penetrate into the magne- 
tized region by a distance ~ 2i?£. cr , while the injected 
electrons (necessary to compensate the CR charge) stay 
close to the injection strip 6 . The plots c) represent the 

6 The fact that CRs can only penetrate by a single distance given 
by ~ 2Rl ct is caused by our particular choice of a monoenergetic 
CR energy distribution. In reality, CRs should have a wide distri- 
bution of energies, so they should be able to penetrate by different 
distances, depending on their energy. We will get back to this point 
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Fig. 6. — Same as Figure|5] but at t-ypcDl = 40. 

CR mean velocities along x and z (pointing out of the 
simulation plane). We can see how the mean CR veloc- 
ity is dominated by its z component (in the three Fig- 
ures of the order of ~ c/2), due to their deflection in 
the By field, which confirms the model presented in the 
lower plot of Figure [I] Finally the plots d) show the 
mean and maximum magnetic amplification as a func- 
tion of x. We can see that the typical wavelengths 
of the loops are consistent with the analytical theory, 
which predicts the dominant wavelength to be such that 
^PCDi/R° L , cr » (vA,o/c) 2 (8irn l /n cr T cr )(5B/B ), where 
we have assumed that the mean CR velocity along z is 
~ c/2. Applying this formula to the cases of Figures 
M |6J and [7] we see that the dominant length scales of 
the fluctuations are about 1.5 times smaller than what 
the t heory predicts. Also, as anticipated at the end of 
§3.2| the charge separation due to loop expansion does 
not stop the non-linear growth of the PCDI in this case. 
The reason is related to the way the background plasma 
compensates the charge of CRs, as they penetrate into 
regions of low initial CR density. We refer the reader to 
Appendix [B] for details on this point. 
A clearer view of the spatial structure of the PCDI can 

at the end of j ]3.3| 
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Fig. 7. — Same as Figure|5] but at t-ypcDl = 44. 

be seen from Figure [8] which shows the components of 
the field and the density fluctuations for a selected region 
of simulation 5* at t^ PCDI = 32. The observed features 
of the PCDI turbulence are essentially the same as for 
the case of the externally imposed current: the non-linear 
field adopts a loop- like structure with very strong density 
contrasts (Sp/po ~ 10). 

As was discussed before, the PCDI happens while the 
background plasma is being pushed by the — J cr x Bq 
force. This pushing can also amplify the field due to a 
simple field compression. In order to differentiate the 
growth due to the compression from the one due to the 
PCDI modes, the dashed lines in Figs. [5] [6] and [7] show 
the results of an analogous simulation where the y di- 
mension of the box has been significantly reduced. This 
simulation behaves as a one-dimensional run, where the 
PCDI modes (whose k\\y) cannot grow. We see that 
the magnetic amplification corresponding to this one- 
dimensional case is significantly smaller than the one due 
to the growth of the PCDI. 

3.3.2. PCDI saturation 

Our simulations show that, similarly to the case of the 
CRCD instability, the field saturation happens when the 
Larmor radii of the CRs, Rl,ct, become comparable to 
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Fig. 8. — The spatial structure of the density fluctuations (a), 

magnetic field strength (6), and the x, y, and z components of B 
(plots c, d, and e) are depicted for a region of simulation S 1 * at 
tlPCDl = 32. The arrows in plot 6) represent the field projectin 
on the plane z = 0. The density and the field are normalized in 
terms of po an( i Boi respectively. 
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Fig. 9. — The left and right panels show magnetic amplification 
(normalized in terms of So) as a function of x at the moment of 
saturation for simulations analogue to S*, but with va o/ c = 1/160 
and 1/40, respectively. The red and blue lines show the maximum 
and average amplification, respectively. And the black line shows 
the magnetic amplitude of a field at local energy equipartition with 
the CRs. 

the size of the generated magnetic fluctuations. In our 
fiducial simulation S*, the magnetic field reaches satu- 
ration at tjpcDi — 40 (Figure pi, where it is amplified 
by a factor of ~ 10. This implies that the Larmor radii 
of the CRs,i?i jCr , should be ~ l/10th of its initial value, 
R° L cr . This coincides with the typical length scale of the 



magnetic fluctuations observed in Figure [6j where the 
distance is measured in units of R° L cr . This saturation 
criterion can be understood in the same way as for the 
case of the CRCD instability. Since the magnetic growth 
is driven by the CR current (in this case perpendicular 

to -Bo) j the growth will stop when this current is signifi- 
cantly reduced, which happens when CR trajectories are 
strongly perturbed by the generated field. This requires 
the length scale of the fluctuations to be comparable to 
Rl,cv We check this in Figure [7J where the z component 
of the CR mean velocity (in regions where n cr > 0.1n cr n) 
has been reduced to almost half of its value in Figure IB] 
By considering the expected evolution of the size of the 
magnetic fluctuations (sa 2ir /kpcDi ~ cB/J cr ) and as- 
suming that the CR mean velocity perpendicular to the 
field is ~ c/2, it is possible to estimate the maximum field 
reachable by the PCDI. If we make cB/J cr Rl,ct = 



rC T cr /e cr B, where 



and e rr are the mass and 



charge of the CRs, its possible to show that saturation 



will happen at B /8n 



r c 2 r cT ./167r. This im- 



plies that the field at saturation should be a factor of 
~ Vl67r 7 below energy equipartition with the CRs. 
We can test this criterion by looking at panels d) of Figs. 
[6jand|7] We see that the maximum amplification is about 
~ 3 times smaller than the one corresponding to local 
energy equipartition between CRs and magnetic field. 
The mean magnetic amplification is an extra factor of 
~ 3 smaller, due to the existence of regions of very low 
density and magnetic field amplitude in the PCDI tur- 
bulence. We can also check this condition by seeing the 
saturation values of the field for simulations with the 
same CR energy density as in simulation S*, but with 
different values of Bo. The left and right panels in Fig- 
ure [9] show the mean and maximum magnetic amplifica- 
tion as a function of x at saturation for simulations with 
v A,o/c — 1/160 and 1/40, respectively. Again, the ob- 
served saturation is consistent with the mean magnetic 
field being ~ 7 times smaller than the one corresponding 
to energy equipartition with the CRs. 
The main uncertainty in this saturation criterion is given 
by the particular CR energy distribution we use, which 
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corresponds to monocnergetic CRs. This implies that 
CRs are able to penetrate into the simulation box by a 
distance shorter than ~ 2i?£. cr . This defines a region 
or "interface" where PCDI modes grow, with the rest of 
the plasma only being compressed due to the accelera- 
tion of the plasma in the interface. The beginning of 
the compression of the central part of the region of pre- 
amplified field can be seen in Figure [7] where the den- 
sity of the background plasma shows prominent spikes 
at x/R° L cr « 2 and 3. However, in reality, CRs should 
have a wide distribution of energies. So, in principle, 
the magnetized CRs (the ones with i?z, jCr < Ao) should 
be able to fill most of the volume of pre-amplified field. 
Also, since the number density of the CRs is expected 
to be a strongly decreasing function of their energy, low 
energy CRs (with small Rl,ct and, therefore, filling a 
small fraction of the volume of the pre-amplified field 
region) will produce PCDI modes on length and time 
scales significantly shorter than the ones of the highest 
energy CRs (that still satisfy the magnetization condi- 
tion, Rl.ct ^o)- Thus, in a realistic setup, the growth 
of PCDI field should start at the boundary of the regions 
of pre-amplified field and then propagate towards their 
interior due to the action of higher energy particles. By 
choosing monoenergetic CRs, we ignore this effect, and 
concentrate on capturing the basic physics of field ampli- 
fication by the PCDI instability, including the growth of 
PCDI modes due to the formation of a J cr perpendicular 
to Bq , and the effect of the CR back-reaction on the final 
saturation of the instability. A more detailed study in- 
cluding the effect of a broad CR energy distribution will 
be presented elsewhere. 

4. APPLICATION TO SNR SHOCKS 

We have found that the saturation of the PCDI is deter- 
mined by the CR back-reaction, and that it would imply 
magnetic fields ~ \/l6n w 7 times below energy equipar- 
tition with the CRs. We can now apply this result to 
the expected conditions in the upstream medium of SNR 
shocks. If we assume that all CR energy is carried by a 
monoenergetic CR population, we can estimate the max- 
imum amplification factor given by 



(6B/B ) r 
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v sh lOkm/s 
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a/s/ V \ 
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(11) 



where rj = u cr /u t h is the ratio between the energies of the 
CRs, u cr , and the one of the downstream thermal plasma, 
u t h- The main uncertainty in this estimate is given by 
the assumption that CRs are monoenergetic, so they will 
only amplify the field in an interface of width ~ 2i?x, lCr 
within the regions ofpre-amplified magnetic field (see 
lower panel in Figure [l I. In reality, CRs will have a wide 
distribution of energies, so they are expected to enter 
most of the volume where the field has been previously 
amplified. Since the density of CRs decreases with their 
energy, the growth of PCDI modes should happen faster 
in the outer part of the magnetized regions, and gradu- 
ally propagate towards the inner parts due to the action 
of the higher energy CRs. The exact details of how this 
process unfolds, and how it would change the estimated 
amplification of Eq. (Ill, requires further study, taking 
into account more realistic CR energy distribution. 



The maximum amplification [TTJ is about 4 times larger 
than for the case of the CRCD instability. Since the satu- 
ration criteria of the two mechanisms are analogous, then 
the factor of ~ 4 is explained by the different mean CR 
velocities used for calculating the corresponding J cr (per- 
pendicular and parallel to Bo for the PCDI and CRCD 
instability, respectively). While for the CRCD instability 
the mean velocity of the CRs can be considered close to 
v a h, i n the case of the PCDI it is ~ c/2. Thus, the max- 
imum amplification should differ by a factor y // (c/2)/v s f l 
(sa 4 for v s h — 10 4 km/s). If the estimate given in Eq. 



(Ill were calculated considering the CRs that are already 
diffusing with a speed of about v s h, then both instabili- 
ties would saturate at about the same level. 
It is important to make sure that the PCDI would have 
time to grow before being advected into the shock. The 
growth of the PCDI field is expected to happen at a dis- 
tance from the shock that is at least the distance that 
CRs protrude into the pre-amplified field, ~ 2i?£ iCr (this 
would correspond to the field loop that is the closest to 
the shock, in Figure fTJ. Then, we can define an ad- 
vection time, t a dv, that would be at least 2R^ cr /v s h- 
In addition to that, the typical time scale of growth is 
given by the inverse of the growth rate of the instabil- 
ity, Ipcdi = {Jcr{ir/ pc 2 )) 1 / 2 . Therefore, the product 

tadvlPCDi ~ 13(77/0.05) (tWl0 4 km/s)/(4«4, /10km/s) 
(assuming CRs moving on average at ~ c/2) shows that, 
as long as the growth of the field starts at a distance a 
few times larger than the typical CR Larmor radius, the 
advection time would not be a restriction to the growth 
of the PCDI. 

The PCDI would also contribute to the efficiency of 
the acceleration of particles at the shock. Indeed, the 
fact that the field saturation happens at Rl,ct ~ ^pcdi 
means that the lower energy CRs would be able to gener- 
ate magnetic fluctuations at length scales comparable to 
their Larmor radii, which would increase their diffusion 
and, therefore, the efficiency of their shock acceleration. 

5. SUMMARY AND CONCLUSIONS 

In this work we proposed a new plasma instability, 
the perpendicular current-driven instability (PCDI), as 
a candidate to amplify magnetic fields due to CRs dif- 
fusing in front of SNR shocks. The PCDI consists of 
purely growing, compressional waves produced by the 
CR current, J cr , perpendicular to the initial magnetic 
field, B . The time and length scale of growth of the 
PCDI are similar to the ones of the CRCD instability, 

i-e., Ipcdi ~ J cr {^l pc 2 ) 1/2 and X PC di ~ cB/J cr . How- 
ever, whereas the CRCD instability is driven by a current 
parallel to B , the PDCI is produced by a current per- 
pendicular to the initial field. Thus, the fastest growing 
instability will be determined by whether J cr is quasi- 
perpendicular or quasi-parallel to the Bq. 
We show that, in the upstream medium of SNR shocks, 
the required perpendicular current can be due to CRs 
with Larmor radii, Rl,ct, smaller than the length scale 
of a previously amplified magnetic turbulence, Ao- This 
scenario is motivated by earlier PIC studies showing that, 
far from the shock, the CRCD instability can produce 
non-linear magnetic amplification, driven by the highest 
energy CRs. This amplification, although limited to a 
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factor of less than ~ 10 for the expected energy carried 
by CRs, would be characterized by non-linear fluctua- 
tions in the plasma density an d magnetic field strength 
( Riquelme fc Spitkovsky|[2009 ) . The typical length scale 
ot these fluctuations would be close to the Larmor ra- 
dius of the highest energy CRs (which is the saturation 
criterion that implies a maximum amplification factor 
of ~ 10), producing a situation where most CRs would 
have a Rl,ct significantly smaller than the typical size 
of the pre- amplified fluctuations. We show that, under 
these conditions, the current perpendicular to B$ can be 
produced by lower energy CRs protruding into the re- 
gions of pre-amplified field by a distance ~ 2i?£ cr , which 
would give rise to a mean velocity of magnitude ~ c/2 
perpendicular to the field. Due to this, in regions close 
to the shock, J cr would be quasi-perpendicular to the 
field, making the PCDI grow faster and on length scales 
smaller than the scales of the CRCD instability. 
We used the two-dimensional PIC simulations to study 
the evolution of the PCDI. We showed that the insta- 
bility exists and confirmed its theoretical growth rate. 
We also studied its possible saturation mechanisms, and 
found that the strong deflection of CRs in the amplified 
field ultimately stops the magnetic growth, which hap- 
pens when the the Larmor radius of the CRs, Rl,ct, is 
close to the size of the magnetic fluctuations, Xpcdi- 
This saturation mechanism is qualitatively the same as 



the one found fo r the CRCD instability ( Riquelme fc 
Spitkovsky||2009 1 . However, since the perpendicular CR 



current is due to a CR mean velocity of magnitude ~ c/2, 
then the maximum magnetic field amplification would be 
~ \A67r ~ 7 times below energy equipartition with the 
CRs. This would imply a maximum amplification fac- 
tor in the upstream medium of the SNR shocks of ~ 45. 
The maximum downstream amplification could increase 
to ~ 200, including the compression at the shock. 
The PCDI would also contribute to the efficiency of the 
acceleration of low energy particles at the shock, by pro- 
viding magnetic fluctuation on scales comparable to their 
Larmor radii, which allows them to diffuse. Given the 
geometry of the PCDI (see Figure [2|, where the fluc- 
tuations in the fields and particle properties happen in 



a plane perpendicular to J cr , we believe that our two- 
dimensional simulations (with J cr pointing out of the 
plane of the simulation) are already capturing the essence 
of the instability. However, future three-dimensional sim- 
ulations will be important to investigate possible extra 
effects that may appear by including the additional di- 
mension. 

One interesting signature of the PCDI is that it predicts 
the amplified field, B, that satisfies B 2 tx v 2 h , assuming 
a weak dependence of r\ (= u cr /itt/j) on w s ^. This rela- 
tionship contrasts with the one ob tained for the CRCD 
instability acting alone, B oc v^ h peUpOOl, and con- 
stitutes a possible observational test to shed light on the 
nature of the magnetic amplification in S NR shocks (see 
observational results of |Volk et al.||2005| . 
Although we have focused on the case of non-relativistic 
shocks in SNRs, the PCDI may also be relevant for 
the case of upstream magnetic amplification in the rela- 
tivistic shocks of jets and Gamma Ray Bursts. In that 
case, as seen from the reference frame of the upstream 
medium, the CRs will only experience a small deflec- 
tion before being advected into the downstream medium 
of the shock. Thus, CRs will provide an electric cur- 
rent parallel to the shock normal that may also amplify 
PCDI modes. The amplification of magnetic field due 
to the expansion of magnetic loops (as in the non-linear 
regime of the PCDI) in the upst ream medium of GRB 
shocks was already discussed by |Milosavljevic fc Nakar| 
(2006). They suggest that the growth of the field would 
saturate when neighboring loops collide. Our simula- 
tions show that these collisions do not stop the growth. 
Instead, they make neighboring loops merge and grow in 
size. We also found that under certain conditions, elec- 
trostatic forces due to the charge separation involved in 
the expansion of the loops may contribute to the quench- 
ing of the nonlinear growth. This possibility requires a 
dedicated study, which will be presented elsewhere. 
In conclusion, we have shown that the PCDI constitutes 
a viable amplification mechanism that, in combination 
with the CRCD instability, strengthens the idea of CRs 
being responsible for the significant fraction of magnetic 
amplification inferred from SNR shock observations. 
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APPENDIX 

APPENDIX A: THREE-FLUID DERIVATION OF PCDI DISPERSION RELATION 

We calculate the dispersion relation corresponding to the growth of magnetic field fluctuations, SB, due to CRs 
propagating at velocity v cr perpendicular to the initial magnetic field, Bq. We focus on the case SB _L Bq and k \\ Bq, 
where k is the wave vector of the Fourier modes of the fluctuations. Ions, electrons, and CR are modeled as three 
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separate fluids with densities, rij, n e , and n cri such that n, + n cr = n e . Electrons are approximated as massless 
particles, whereas ions and CRs have finite masses and sound speeds c St i and c SjCr , respectively. 

The essential part of the calculation consists of determining the perturbations to the plasma currents, 5 J, expected 
in the presence of the amplified field. To first order, the current perturbation contributed by species "j" is given by 
SJj = SrijVj + rijSvj, where 5rij and Svj represent the perturbations to the density and velocity of species "j" . In the 
case of the CRs, these quantities can be calculated from the perturbed part of the Euler equation, 

T^i^cr^cr ^ Cs.cr^li^ cr^ ' 8n cr -(- GTl cr SE , . * 

+en cr ( s -^xB +^x5B), 1 ' 

where nii is the mass of the ions (assumed equal to the mass of the CRs) , r cr is the Lorentz factor of the CRs, e is 
the charge of the ions (also, equal to the one of the CRs), and 8E corresponds the the electric field fluctuations. For 
the kind of fluctuations we are interested in, the relevant currents will be the ones perpendicular to Bq- If we assume 
that Srij, Svj, 5E, and 5B are stationary waves that grow exponentially at a rate 7, then 

„ e (78E + u CtCr 5E X y) 

Sv cr , ± = — — ^ 5- — '-, (A2 

where the subscript _L stands for "perpendicular" , ui c . cr is the cyclotron frequency of the CRs, and y is the axis parallel 
to B . 

The CR density can be obtained from the y component of Equation (Al I. If we neglect the term proprotional to 5E, 
we obtain 

—en cr d y (— x 5B) y 
SUcr = mTW + fc 2 ^)' ( } 

This approximation is justified if we look at the y component of the Euler equation for electrons, 

m e n e ^^ = -en e {6E y + -j x SB). (A4) 

Since they are assumed to have zero mass, then SE y — — x 8B. The unperturbed electron velocity v e is given mainly 

by the E x Be/ B 2 drift of the background plasma, which we make initially zero by choosing the appropriate frame of 
reference. 

In the case of ions and electrons, the velocity fluctuations are calculated from the fluctuations of their drift velocities. 
These drifts are the E x Bc/B 2 drift, which we call and an extra drift due t o the time variation of the fields, 
= —(rrijc/ejB 2 )(dvo/dt) x B (see Appendix A of Riquelme fc Spitkovsky||2009 1 . Thus, it is possible to show that 



^0, = ^B Cj (A5) 

and 

(11 rrijC 2 d8E , . . 

where bvf^ affects equally the ions and the electrons, while bv^ 1 only affects the ions (given that electrons are massless). 
The density fluctuations of the ions can be obtained directly from the y component of their Euler equation: 

«M7 + k c Syl ) oy 

It can be shown that in the regime in which |{^°'| v cr , the time variation of the zeroth order component of the 
electric field is given by dE/dt = —Airu cr n cr ev 2 A /c 2 , where va,o is the initial Alfven velocity of the background plasma. 
Thus, the homogeneous part of can be expressed in terms of the time variations of the electric field to obtain, 

x en cr d(% x SB) y 

^= mi(72 + fc2<i ) By (A8) 

Since we are choosing a frame where 6E y is initially zero, the term d y 8E y can also be neglected. Thus, from the 
Poisson equation we obtain that Sn e 1=3 Jrij + 5n cr . This way we have all the density and velocity fluctuations needed 
to calculate 5 J. If we are interested in the case where the CR back-reaction can be neglected, the value of 8 J can be 



Perpendicular Current-Driven Instability 



13 



0-000 
0.9 





7 ' 






) J ) 


V 




A * ' JgF" • 






/ V 





E/3c 



0.267 





Fig. 10. — The spatial structure of the electric field produced by the PCDI is depicted for a region of simulation S» at t^/pcDl = 32. 
The electric field is normalized in terms of Bq. 



calculated by taking the limit T c 



oo. In that case we obtain that 

2 „ m 2 

<■ ^ 

V 



X T 1 C ' St ! en cr a (V, 

SJ= - T j6E —d y - 

47T v A rrtirii V i 



x SB) , ., , , , 

r / y V7 ^ + fe^cf , 



(A9) 



Thus, combining the Ampere's and Faraday's laws with Eq. (A9l we recover Eq. d8|. 



APPENDIX B: THE EFFECT OF CR CHARGE COMPENSATION ON THE GROWTH OF THE PCDI 

If CRs are made of only positively charged particles, then the background plasma is expected to carry a negative, 
compensating charge that nearly cancels it. In j ]3.2| we saw that this background charge may suppress the non-linear 
growth of the PCDI. In regions where the PCDI produces a decrease in the plasma density, a net positive charge 
will grow. In the linear regime, this net charge growth is not important. It creates an electric field proportional to 
the product of the charge fluctuations in the background plasma, which is approximately equal to p c ,crSp/ p, times 
the characteristic length scale of the fluctuations, which is proportional to 5B/B . Thus, the force on the plasma 
due to this effect is of second order in perturbation theory, and can be neglected. In the non-linear regime, however, 
this electrostatic force can be relevant. As the PCDI magnetic loops expands, the electric field in the loops will 
be given by E ~ Ai:n cr ef, where r points in the radial direction. Thus, the loops will feel an electrostatic force 
~ — Aitrn 2 cr e 2 , that will push them back to their center. Recalling that, when the amplified field becomes non-linear, 

\f(t) | ~ Xpcdi ~ cB/J cr , then this force will be of the order of J cr x B if the CRs move at an average velocity close to 

c. This implies that the force driving the expansion (— J cr x B/c) would be roughly compensated by this electrostatic 
force, which would quench the non-linear growth of the PCDI. 

This situation, however, would not happen when t he P CDI is driven by CRs protruding into regions of pre-amplified 
field, as suggested by the simulations presented in S3.3 The reason has to do with the way the compensating charge 
is distributed in the loops. As CRs penetrate into the pre-amplified field regions, their electric current induces a 
return current in the plasma (J re t ~ Jcr), which is what compensates the charge. However, if this happens while the 
loops are being formed, this return current will tend to separate charges within the loops, producing a dipolar charge 
distribution in them. 

This can be seen from the shape of the electric field, E. shown in Figure 10 which depicts the same simulation 
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region as Figure |8| In the magnetic loops locations, E does not have the monopolar shape needed to stop the loop 
expansion. Instead, we see that, whereas the y component has a monopolar shape, the x component looks li ke a 
dipolar field, indicating the presence of charge separation within the loops. In the situation depicted in Figure |10[ 
electrons accumulate to the right of the loops. Thus, while the negative charge of the plasma concentrates mainly in 
regions where B\\y, the places where B _L y will be more neutral. Thus, there will be parts of the loops where the 
electrostatic forces are less importan t th an in the cases where the charge in the loops is distributed homogeneously 
(like in the simulations presented in { 3.2 1, and their evolution will be dominated by the — J cr x B force that expands 
the loops. This way, the PCDI can grow to non- linear amplitudes, even when the CRs may be mainly composed of 
positively charged particles. 



